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Baryon Stopping as a new Probe of Geometric Scaling 
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We suggest to use net-baryon rapidity distributions in central relativistic heavy-ion collisions 
at SPS, RHIC and LHC energies in order to probe saturation physics. Within the color glass 
condensate framework based on small-coupling QCD, net-baryon rapidity distributions are shown 
to exhibit geometric scaling. In a comparison with RHIC data in Au + Au collisions at ^/snn = 62.4 
GeV and 200 GeV the gradual approach to the gluon saturation regime is investigated. Predictions 
for net-baryon rapidity spectra and the mean rapidity loss in central Pb + Pb collisions at LHC 
energies of ^/sjvjv = 5.52 TeV are made. 
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Baryon stopping in relativistic heavy-ion collisions as 
a probe of QCD-matter at high parton density is of great 
current interest [H, H, El, Hj] . Theoretical QCD-based ap- 
proaches usually focus on charged-hadron production. In 
the central rapidity region a reasonable understanding 
has been achieved in the color glass condensate (CGC) 
framework 0, S S @] through inclusive gluon production 
@; Ell • In this theory, due to the self- interaction of glu- 
ons, the number of gluons in the nuclear wave function 
increases with increasing energy and decreasing longitu- 
dinal momentum fraction x carried by the parton. 

Unitarity requires that the gluon density saturates be- 
low a characteristic momentum scale, the so-called sat- 
uration scale Q s . In this regime gluons form a coher- 
ent state. Presently the evidence for the existence of 
this state of matter is, however, not yet clear. Due to 
the dependence of the saturation scale on rapidity and 
mass number, it has been proposed that saturation effects 
should be studied with heavy nuclei and large rapidities 
at RHIC energies and beyond. 

In this Letter we suggest to use the rapidity distribu- 
tion of net baryons (B — B) in central heavy-ion colli- 
sions as a testing ground for saturation physics, cf. Fig. 
[TJ In A + A collisions, two distinct and symmetric peaks 
with respect to rapidity y occur at SPS energies [111 ] and 
beyond. The rapidity separation between the peaks in- 
creases with energy, and decreases with increasing mass 
number A reflecting larger baryon stopping for heavier 
nuclei, as has been investigated phenomenologically in 
the relativistic diffusion model [121 ]. 

The net-baryon number is essentially transported by 
valence quarks that probe the saturation regime in the 
target by multiple scatterings [HI]. During the collision 
the fast valence quarks in one nucleus scatter in the other 
nucleus by exchanging soft gluons, leading to their redis- 
tribution in rapidity space. We take advantage of the fact 
that the valence quark parton distribution is well known 
at large x, which corresponds to the forward and back- 
ward rapidity region, to access the gluon distribution at 
small x in the target nucleus. Therefore, this picture pro- 
vides a clean probe of the unintegrated gluon distribution 
ip(x,pr) at small x in the saturation regime. Here pt is 



the transverse momentum transfer. 

We have two symmetric contributions, coming from 
the two beams. The contribution of the fragmentation of 
the valence quarks in the forward moving nucleus is given 
by the simple formula [l4| for the rapidity distribution of 
hadrons: 



dN 
dy 



C 



d 2 Pr f n , , \ 
— — xiq v (xi,Qf) ip(x 2 iPT) 

Pt 



(1) 



where xi — pt / t/s exp(y) , X2 — pt / \/s exp(—y) are the 
longitudinal momentum fractions carried, respectively, 
by the valence quark in the projectile and the soft gluon 
in the target. The factorization scale is set equal to the 
transverse momentum, Q f = pt- The contribution of 
valence quarks in the other beam nucleus is added in- 
coherently by changing y — > —y. The gluon distribu- 
tion is related to the forward dipole scattering amplitude 
Af(x, 7"t), for a quark dipole of transverse size tt, through 
the Fourier transform 



<p{x,Pt) 



2ttp t 



r T dr T J\f(x,r T )Jo(rTPT)- (2) 



In the fragmentation region of the projectile the valence 
quark parton distribution function (PDF) is dominated 
by large values of x\. We integrate out the fragmenta- 
tion function such that the hadron rapidity distribution is 
proportional to the parton distribution. The overall con- 
stant C depends on the nature of the produced hadron. 

One important prediction of the color glass conden- 
sate theory is geometric scaling: the gluon distribution 
depends on x and px only through the scaling vari- 
able p T /Q 2 s {x), where Q 2 g (x) = A 1 ^Ql x ~ x , A is the 
mass number and Qo sets the dimension. This has been 
confirmed experimentally at HERA [15]. The fit value 
A = 0.2 — 0.3 agrees with theoretical estimates based 
on next- to- leading order Balitskii-Fadin-Kuraev-Lipatov 
(BFKL) results [H,[i3]- To show that the net-baryon dis- 
tribution reflects the geometric scaling of the gluon dis- 
tribution, we perform the following change of variables: 
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FIG. 1: (color online). Rapidity distribution of net baryons 
in central (0 - 5%) Pb + Pb collisions at SPS energies of 
y^ijvjv = 17.2 GeV (top frame). The theoretical calculations 
are compared with NA49 results that have been extrapolated 
from the net-proton data Dashed curves are for A = 

and Qq = 0.08 GeV 2 , solid curves are for A = 0.15 and Q 2 — 
0.07 GeV 2 , and dotted curves are for A = 0.3 and Q 2 , = 0.06 
GeV 2 . At RHIC energies of v /57v7v = 62.4 GeV (middle frame, 
- 10%) and 200 GeV for central Au + Au, our corresponding 
theoretical results are compared with BRAHMS net baryon 
data (circles) [2, 0]. At 200 GeV, triangles are preliminary 
scaled BRAHMS net proton data points for - 10% PH. 



Thus, we rewrite Eq. ([I]) as 

^W = f T-^W^+V), (4) 

dy 2tt J x 

where r = Ii^s/Qq) — In A 1 / 3 — 2(1 + A) y is the corre- 
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FIG. 2: (color online). Rapidity distribution of net baryons in 
central Pb + Pb collisions at LHC energies of ^/sjvjv = 5.52 
TeV. Theoretical distributions are shown for A = (dashed), 
A = 0.15 (solid), and A = 0.3 (dotted curve), with Qq - values 
as in Fig.l. 



sponding scaling variable. Hence, the net-baryon multi- 
plicity in the peak region is only a function of a single 
scaling variable t, which relates the energy dependence 
to the rapidity and mass number dependence. In the 
fragmentation region, the valence quark distribution is 
only very weakly dependent on Q/. 

From the equation for the isolines, r = const, one gets 
the evolution of the position of the fragmentation peak 
in the forward region with respect to the variables of the 
problem 



1 



2/peak — 



1 + A 



Vbc 



In A 1/6 ) + const, (5) 



where y bc am = 1/2 • ln[(£ + pl)/{E - p L )] ^ hxy/sfm 
is the beam rapidity at beam energy E and longitudinal 
momentum pl with the nucleon mass m. 

To take into account saturation effects in the target 
we choose the Golec-Biernat-Wiisthoff model [l9| for the 
forward dipole scattering amplitude Af, leading to (cf. 
Eq. © and H3) 



(p(x,p T ) = Air 



Pt 



exp 



o 

Pt 



(0) 



The valence quark parton distribution of the nucleus is 
taken to be equal to the valence quark PDF in a nucleon 
times the number of participants in the nucleus. We are 
focusing here on the forward rapidity region, and interpo- 
late to mid-rapidity where small- a; quarks are dominant, 
by matching the leading-order distributions [2(| and the 
Regge trajectory, xq v oc x - 5 , at x = 0.01 

To account for large-x effects in the gluon distribu- 
tion, we multiply the distribution function by (1 — X2) 4 
[l(| • Mass effects are considered through the replacement 

p T -> y / p T + m 2 - 

Our results for net-baryon rapidity distributions in 
central Pb + Pb and Au + Au collisions are shown in 
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Fig. Q] Dashed curves are for A = 0, solid curves for 
A = 0.15, and dotted curves for A = 0.3, with the corre- 
sponding Qp-values fixed at SPS, as given in the caption. 

We compare with SPS NA49 Pb + Pb results at 
y/sW^ = 17.2 GeV [U, and BRAHMS Au + Au data at 
62.4 GeV and 200 GeV [EH El]. We obtain the number 
of baryon participants at SPS energy in the full rapid- 
ity range from a double-gaussian fit of the NA49 data 
for - 5% central Pb + Pb collisions as N B = 380. The 
normalization in our model calculation is 12% lower than 
this value since we do not account for the baryons near 
2/bcam in the tails. We maintain this correction at RHIC 
energies where the tails are in the unmeasured region. 

For Au + Au at RHIC energies, we take Glauber 
results for the number of participants: At 62.4 GeV 
N B = 314 ± 8 for centrality - 10% @ , and at 200 
GeV N B = 357 ± 8 for - 5% fl. The comparison with 
the SPS and RHIC data slightly favors A < 0.15 and 
hence, the asymptotic regime with A ~ 0.3 is not yet 
reached at RHIC. 

Our prediction for central Pb + Pb at 5.52 TeV LHC 
energies is shown in Fig. [2] for A = 0,0.15, and 0.3. At 
LHC energies the mid-rapidity region is almost baryon 
free, we obtain dN/dy(y — 0) ~ 4 for net baryons. The 
position of the fragmentation peak is very sensitive to 
the value of A, with a difference of about 1.5 units of 
rapidity between the A = and 0.3 cases. It is possible 
that the full scaling regime with A approaching 0.3 can be 
reached at or beyond LHC energies, but presently none of 
the LHC-experiments is capable of measuring identified 
protons or neutrons from central Pb + Pb collisions in 
the region of the fragmentation peaks. This would be 
a relevant proposal for future extensions of the detector 
capabilities at LHC. 

Physically the two peaks represent the result of the 
scattering of the fast moving projectile valence quarks in 
the target, they are deflected, their distribution broadens 
and carries information about the gluon distribution in 
the target. This is in analogy to x-rays that are deflected 
by a crystal and carry information about its structure. 

With increasing energy the peaks move apart, the solu- 
tions behave like travelling waves in rapidity space [2l| . 
which can be probed experimentally at distinct values 
of the beam energy, or the corresponding beam rapidity. 
We have derived the peak position as a function of the 
beam rapidity as j/ pca k = v j/bcam + const with the peak 
velocity v = 1/(1 + A), cf. Eq. (JSJ. The position of 
the peak in rapidity space as a function of the beam en- 
ergy can in principle be determined experimentally, or at 
least estimated (RHIC). Theoretically, its evolution with 
energy provides a measure of the saturation scale expo- 
nent A. Hence, a precise determination of the net-proton 
fragmentation peak position as a function of beam en- 
ergy would provide detailed information about the gluon 
saturation scale. 

In Fig. [3J we show our numerical results for the mean 
rapidity loss (Sy) = j/bcam — (y)- At low energies they 
agree with the experimental values of baryon stopping 
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FIG. 3: (color online). The mean rapidity loss (Sy) as ob- 
tained from our theoretical results is plotted as a function of 
beam rapidity {/beam, solid curve. The star at j/beam = 8.68 is 
our prediction for central Pb + Pb at LHC-energies of y'sjvjv 
= 5.52 TeV with A = 0.15, the dashed curve is for A = 0, 
the dotted curve is for A = 0.3, with Qq values as in Fig.l. 
Analysis results from AGS Au + Au data (E917, E802/E866, 
triangles) [H, SPS Pb + Pb data (NA49, square) [U, RHIC 
Au + Au data (BRAHMS, dots, with triangles as lower and 
upper limits) [l], Q] are compared with the calculations. 

that have been obtained at AGS and SPS energies pH[22| 
irrespective of the value of A. Here we have considered 
the effect of the missing particles in the tails as described 
above for the rapidity distributions, placing them halfway 
between the mean rapidity, and the beam value. At 
RHIC energies of 62.4 GeV and 200 GeV, the mean ra- 
pidity loss depends on A, and the A = 0.3 result (dot- 
ted curve) is beyond the upper limit given by BRAHMS 
[HQ, whereas A < 0.15 is consistent with the upper limit 
of the data. Consequently, up to the highest RHIC ener- 
gies the expected scaling regime with A ~ 0.3 [15j is not 
yet fully reached, in accordance with [Tol. [24l|. 

Our result emphasizes the importance of a detailed 
measurement at LHC energies to allow more definite con- 
clusions about the value of A, which would then be de- 
termined by the slope of the mean rapidity loss at high 
beam rapidity above RHIC (solid curve in Fig. 3). 

Assuming that the mean rapidity evolves similarly to 
the peak position, (y) = y pea k-l-const., the linear increase 
of the mean rapidity loss at large energies corresponding 
to beam rapidities ybcam > 5 is given by 

^) = -i ^ x i/bcam + COnSt. (7) 
t + A 

Hence, the mean rapidity loss that accompanies the en- 
ergy loss in the course of the slow-down of baryons pro- 
vides at large beam rapidities ybeam > 5 a measure for A 
and thus, a test of saturation physics. The case A = 0, 
or equivalently Q s constant, leads to a saturation of the 
mean rapidity loss at high energies, and correspondingly 
at large beam rapidities. 
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In the peak region, the average x in the projectile is 
x ~ 0.2 — 0.3, which corresponds to the average momen- 
tum fraction carried by a valence quark. In the target, 
x = (0.2 — 0.3) e~ 2j,poak , it decreases with increasing en- 
ergy. In this kinematic regime we have a natural intrinsic 
hard momentum, the saturation scale Q s - This justifies 
the use of small-coupling_techniques in QCD for calcu- 
lating integrated yields [23j . The effects of the medium 
are expected to be small at forward rapidity since the 
fast moving valence quarks escape the interaction zone 
quickly. A detailed measurement of the peak region 
would then enable us to reconstruct the gluon distribu- 
tion from Eq. {J). 

To summarize, we have presented a saturation model 
for net-baryon distributions to investigate the gradual ap- 
proach to the gluon saturation regime at RHIC energies 
and beyond. In a comparison with BRAHMS net-baryon 
results for central Au + Au collisions at 62.4 and 200 GeV 
we have determined a saturation-scale exponent A < 0.15 



and hence, the full scaling regime is not yet reached at 
RHIC. This result is in agreement with studies of particle 
production that point out a slower growth of the satura- 
tion scale at RHIC energies than the HERA estimate of 
A ~ 0.3 suggests (MSI- 

In particular, we have shown that the peak position 
in net-proton rapidity distributions of centrally colliding 
heavy ions at ultra-relativistic energies obeys a scaling 
law involving the atomic mass and the beam energy. Our 
result for the mean rapidity loss in x/sjvjv = 62.4 GeV 
and 200 GeV Au + Au is for A < 0.15 consistent with the 
upper limit of the corresponding BRAHMS experiments. 
We emphasize the importance of a detailed analysis at 
LHC energies. 
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